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Vector potential of the viscosity wave 
Fiber radius 

Adiabatic speed of sound 
Fiber separation (in the model of a 
regular absorber) 

Amplitudes of wave types in 
air 

Amplitude of temperature wave inside 
the fiber 

Internal wave conductance of the absorber 
Hankel functions of n-th order, 
second type 

Porosity of the absorber 
Fill factor 

Wave numbers of wave types 
Fiber density (number of fibers per 
surface uni* ) 

Stationary pressure fraction 
Alternating pressure 
Radius of an elementary cell (fiber 
in the tube) 

Cylinder coordinates 
See equation (151) /sic/ 

Stationary temperature of the air 
Alternating temperature of air 
Stationary temperature of the fiber 
Alternating temperature of the fiber 
Speed of sound 

Internal wave resistance of the 
absorber 

Surface impedance (with #: standardized 
to (111 co)- 

General cylinder function of 
n-th order 

Temperature conductance coefficient of air 
Temperature conductance coefficient of the 
fiber material 

Propagation constant (with #: standardized 
to '"/'■«), ) 

Adiabatic exponent 
Dynamic viscosity 
Scalar potential 
Scalar potential of 
Scalar potential of 
Thermal conductance 
Thermal conductance 
Kinematic viscosity 
Angular frequency 
Stationary density of air 
Alternating density of air 
Specific flow resistance 
Flow resistance 


the compressional wave 
the temperature wave 
of air 

of the fibre material 


Indices : 


a Temperature wave in air 

i Temperature wave of the fiber 

v Viscosity wave 

i> Compress ional wave 

i (prime): emitting field 

" (double prime): scattering field, constant field 
- (circumflex): amplitudes determined by model geometry 


A MODEL THEORY FOR THE FIBROUS ABSORBER, 

PART 1: REGULAR FIBRE ARRANGEMENTS 

F. P. Mechel 

University of the Saarland Acoustics Laboratory, Saarbrucxen 
1. Introduction 
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Most sound absorbers are fiber absorbers: the absorption materi- 

als consist of loosely interconnected fibers. The fibres themselves 
are usually glass or mineral fibers, but also organic fibers and — for 
special applications — metal fibers are used. 

To describe sound propagation in these fiber materials, at present 
there are practically two theories. The first, which we will call the 
theory of quasi-homogeneous absorber, describes the absorber as a homo- 
geneous, isotropic medium and considers the losses in sound energy in the 
absorber by a flow resistance and tries to describe the structure by a 
"structure factor" in the force equation. The second theory replaces 
the absorber by a model of a bundle of parallel tubes (Rayleigh model) 
with sound resistant, infinite heat conducting walls. Both theories 
have their obvious and known deficiencies. 

In an earlier study /!/ it was found which of the two theories 
was best suited for further development. It turned out that the theory 
of the quasi-homogeneous absorber leads to internal contradictions in 
addition to the disadvantage of frequency dependent material constants 
which cannot be eliminated by a simple change in the theory. Conversely, 
the theory of the tube model appears to correctly reproduce the physical 
processes of sound propagation in the fiber absorber, but suffers by 
having excessive deviation between model and structure and real absorber 
structure . 

Therefore it was suggested to apply the principles of the theory 
of the Rayleigh model, namely the solution of the most exact and com- 
plete differential equations to a model better adapted to the real 
absorber structure. 


*Numbers in the margin indicate pagination in the foreign text. 


The model used here consists of parallel fibers of the same 
material just like the real absorber. The solutions of the complete 
differential equations (L.e. under consideration of air viscosity be- 
tween the fibers and thermal conductance of both the air and fibers) 
are given. The sound wave should propagate In the direction of the 
fibers . 

Initially omitted properties of a real abosrber are: 

1) scattering of fiber radii about an average value 

2) melt beads and adhesive globs 

3) statistic orientation of the fibers In a plane. 

While work with this theory was underway, the excellent papers 
by Attenborough et al. /2 /, fU appeared. There, the theory of multiple 
scattering from cylinders was applied to the fiber absorber. This Is 
a method tailored specifically to the absorber attacked perpendicular 
to the fiber direction. 

Since we are assuming a sound propagation parallel to the fiber 
orientation, the model described here represents a supplement to the 
work of Attenborough. In addition, the calculation presented here 
represents a logical continuance of the Rayleigh model. In the Rayleigh, 
model the sound propagates axis-parallel In the tube. We eliminate the 
main error of this model, namely the tube, whose wall consists of fibre 
material, and study Instead the axis-parallel acoustic propagation In 
a fiber bundle. 

In an Initial equation for this model we assume a regular arrange- 
ment of fibers In a quadratic grid. This is certainly a flaw In the 
model by comparison to a real absorber. Its effect Is even more omi- 
nous since we frequently had to use symmetry properties of this ar- 
rangement In the course of the calculation. The method of calculation 
shows which changes occur In the result If Instead of quadratic symmetry, 
a symmetry of higher order is postulated. 


On the basis of these results a simplified model is discussed 
in a subsequent part of this work (part II). Here, each fiber is 
surrounded by a imaginary cylindrical symmetry surface, i.e. we use 
a symmetry of infinite order. As with the Rayleigh model, we only 
need to find the solutions in this type of cylinder with one fiber 
on the axis. 


2 . Basic Equations /55 

The scalar quantities of state of the sound field, namely, pres- 
sure p , density t> , temperature T are split apart into stationary 
fractions with the index zero and into chronologically variable frac- 
tions having the index one in accordance with exp (j«t). .For the 
sound field quantities the usual linearization assumptions are made 
that their squares are negligible compared to linear terms. We then 
have : 

the equation of force (Navier-Stokes ) : 

^ I 

( V -1- mwl in — i’ An — 
ft oo 

— J rf'i'iul'liv v ----- 0 , 

the continuity equation: 

<\<1 

(), ( 2 ) 

the equation of thermal energy: 
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the equation of state: 
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and the equation of thermal conductance in the fiber: 
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In equation 1 we neglect convective acceleration and in equation 
3, thermal transport by radiation. 


3. Potential Functions and Wave Equations 


The standard potential functions for the acoustic field are de- 
fined by : 


" ^ni(l0 | rot A (d) 

with the secondard condition 

< I i v A 0. (■) 


After substituting into equation 1 under consideration of the time 
law exp (jut) and the identity: 

A f'rmltliv -rot. rot 

rot = red 


we obtain: 
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Since both summands are independent of each other, they must disappear 
individually and we have: 
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Equation 10 comprises the first wave equation, namely for the viscosity 
wave : 




with : 
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Substituting equation 6 Into equation 2 gives: 

(M (12) 

l/o j 

and substituting equation 6 Into equation 3 gives: 

jo; Y’i - - (y ■ I) 7'<> A0 jiA'/'i 0. (13) 

In order to obtain a wave equation for 0 , we must eliminate the middle 

term in equation 9. Substituting equation 12 into equation 4 gives: 
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Substituting equation 14 into equation 9 gives: 
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Substituting A'/; from equation 14 into equation 9: 
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i-A-0 0 . 
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Equation 16 minus equation 15 multiplied by yu/a results in: 
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with 
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According to Lord Rayleigh, equation 17 can be written as a product 
of two wave equations 

(A | //;)( A | ii) ( l> -0. (l!)) 

z 2 

Here, k fe , k£ are the negative solutions of the quadratic equation by 
A, obtained by equating the brackets in equation 17 with zero. Prom 
this we have: 
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Equation 19 is solved if we set: 


'/> % I , 

(A I ^) f /v o, 

(A | /’;)</>* o. 


( 21 ) 

( 22 ) 
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If we consider the order of magnitude of the constants, we see 
that by approximation we only need to consider the first term in 
parenthesis in front of 4 and 4* in equation 17. Prom this we have: 
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and if we extract the root : 




(">A‘ »)“, /’£ i=t jy«)/a. 


( 21 ) 


K 0 is the wave number of the standarized sound wave in air, that is, 
a compressional wave; K a » is the wave number of a temperature wave. 


Accordingly, three wave types appear in the absorber: the normal 
compression wave defined by ^ , a temperature wave described by 0 a and 
a viscosity wave with vector potential 4 . 


4 . Quantities of State 


It is our objective to find solutions for the potential functions. 
In order to describe the acoustic field the acoustic field quantities 
must be expressed by these potential functions. 

The speed of sound is obtained from equation 6 and 21: 


v ■■ - «rml(f/y I - '/>,,) | mt./f . (2. r ») 

The relative change in density is obtained from equation 12 with equation 
22 and 23 as: 

2 ~ L L ^ 0t ’ 1 ^ 0aJ - ( “ (i) 

Prom equation 9 we obtain the relative sound pressure as: 

* • I If, |- // 7 0 a (27) 
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with the coefficents: 
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As in equation 28, expressions with the indices will appear below. We 
intend to. use these as abbreviations for two equations, each of which 
would otherwise have only one index. Therefore, equation 28 stands 
for: 
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Prom equation 4 we then obtain the relative temperature change as: 
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with the coefficients : 
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The following discussion where the explicit dependence on the kinematic 
viscosity i/ is eliminated also pertains to the coefficients //(.,* and : 
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5 . Boundary Conditions 


/57 


The emission conditions and certain symmetry conditions, depending 
on the model, must be met. The following conditions must be met at the 
fiber surfaces by the resulting mathmatical descriptions of the poten- 
tials: disappearance of the speed with all components, equality of 
temperatures of air and fibers, equality of radial thermal flux density 
in air and fiber at the fiber surface, i.e.: 


t>r , Vq, , V z -2 o 

Ti ~ T U 

. a a 
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As a simplified thermal boundary condition we studied the disappearanc 
of alternating temperature at the fiber surface, l.e.: 

V, (/• a) 0 , 

Although the boundary conditions to be met by the Individual fibers 
are always the same, with the individual model there result various 
descriptions of the potentials (e.g. due to different symmetry pro- 
perties) which are to meet the above boundary conditions. 

6. Individual Fiber 


We first consider an isolated fiber. Let it be embedded in a 
cylindrical coordinate system (r, ( ;> ,z), with the z-axis coinciding 
with the fiber axis. Let us call the fiber radius a. As the time 
function we select exp (j wt), as the z-function exp (— r'z). From 
the geometry of the arrangement, we have immediate axial symmetry, 
i.e. <V<V o, From equations 22 and 23 we have: 

/ a* i a „ \ 

V<v *•' ,• ,-v H) - W 

If we set : 

A = (Ar,A v ,A t ), (37) 


Then because of 


we have : 
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to equation 38 and 
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to equation 39 gives 
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.1 r 0, 


And with equation 11 there results: 

o, 

Prom equation 38 and 39 we have: 


and 

that is 



o 



o, 


'h o. 


Therefore, of the vector potential there remains only the component A (/ , 
for which the wave equation applies: 


/ 08 l 0 
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As solutions of equation 36 or 40 , Hankel functions of zero or first 
order are used. Because of the emission conditions only the second 
type of Hankel function come into consideration if we choose the imaginary 
fraction of the cross-propagation constant to be negative: 


('•)•:<>• (.ii) 

Thus we have : 
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With these equations the boundary conditions are met at the fiber sur- 
face . 

As a result of axial symmetry, ■<.„,(/•.-«) o is also true. Prom 

the other conditions there result the following equations: 

v r (r . - a) 0: 

^ /''V< I'c'H i 0v«) -|- Hr I'll i (/■>«) 0. (-Ml) 


The sum !i'ln equation 46 Is an abbreviation for the sum of two terms. 

CV* 

The first summand Is the one behind the summation sign, the second Is 
obtained from the first If we replace the Index p by the Index a. 

Prom the condition v, (r * z) = 0 we obtain: /58 

| A’, //o(/>") 0. (17) 

I'** 

Prom the equality of the internal and external temperature at the 
fiber surface we obtain: 

'ZMelloM -FiMn") 0 ( 4K ) 

II,* 

and from the equality of the radial heat flux it follows : 

A’| (fj«) o. (4!t) 

V.l '* 

the following expression for the wave equation: 

(A | k'f) T\\ o (BO) 

— jw/ai (> r >l ) 

equation : 

Tti(r.z) 

i "Mnr) m 

I 0 

if <7 I !'-■ (p») 

To determine the propagation constants r the determinant of the equa- 
tion system (46) to (49) is eliminated for the amplitudes E p s E a , E. l/S P } . 
This system of equations can now be written for the sake of clarity as 
follows : 


Prom equation 5 we have 


'with : 


We obtain the solution 


with : 


^ | A,’,. / ' // 1 (/•■,, n ) o, 
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0 V r v II x (r v a) A’j i’\>I 0 , 
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The numerical equation shows that for fiber radii as used In real 
absorbers, the propagation constant of an Individual fiber Is not 
affected. But In order to obtain calcuable deviations In propagation 
constant from the free field wave number, fiber radii of at least 
20 cm must be assumed. 

The calculation performed here permits us In future models to 
separate the field of an Individual fiber from the total field In the 
absorber. 

7 . Quadratic Arrangement of the Fibers 
7.1 Geometry 

The geometric conditions are illustrated In Figure 1. We call 
the fiber radius a, the distance of the axes of two neighboring fibers 
Is d. 


In the cartesian coordinate system(x, y) the fiber axes have 

coordinates (md, nd), whereby m, n are whole numbers. The individual 

fibers are identified by the indices m, n. In the cylinder coordinate 

system (r, <i ), the fiber (m, n) has coodinate (R ^ ). 

mn y mn 


Figure 1: 



point 


P. 


The potentials are periodic in y by njl : 

<l>(l\ r ,z) -.0 I-tt/2,3). 


;»n 


The potentials are symmetric at , r 0 , 71 / 4 . 

In the treatment of this model we assume that the field about a 
fibre is composed of a fraction as would be generated by the Individual 
fiber (emitted field) and a fraction composed of the sum of the emitted 
fields of all other fibers (scattering field). The emitting field is 
assumed to be axial symmetric, the scattering field is dependent on 
the angle. 

Prom the shape of the emitting part of a fiber (see equation 42 to 44) 
there results that this field is inconstant on the axis of this fiber. 

The scattering field caused by the other fibers is naturally constant 
and therefore instead of the concept "emitting field" and "scattering 
field" we also use the expressions "inconstant field" and "constant 
field" . 


7.2 Determination of the Scalar Potential 


For the scalar potential we use the wave equation: 


As solution set we select: 


i a / r \ i fa fa | 

r <> \ <> / ' l " /■“ y 1 <V“ 1 x 
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This equation already considers the symmetry about y <». 
Because of the periodicity all terms where i is not divisible by 
four are eliminated in the sum, i.e.: 

CO 

I £ 'C. I : J u U'v.a.r) ci >s (4 / y )| . (55) 
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Here, the amplitudes of the emitting field have a single prime, those 
of the scattering field a double prime. Equation 55 meets the wave 
equation 


i fi / ? \ /a 


<> 

( 50 ) 


with 


4* I 


< o. 

It also meets the conditions of axial symmetric excitation of the 
emitting field, emission condition, the periodicity properties of 
the scattering field. 

Now the scattering field should be described as the sum of the 
emitting fields of other fibers, i.e., we should have: 

oo 

2 /C -hi (<■>•) cos (4 i</) (r>7) 

I (I 

-- h ^ (l % | I'wn /' 0OS(yi„i H </'))• 


Here we have : 

K mu d\ m* !■»», 

cos t/j„j ,i — w/|' w- | //- , (^H) 

siiw/i,„ „ w/l'wi® |- . 

according to the Gegenbauer addition theorem, the right side of equa- 
tion 57 is: 
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In equation 59 we know: 


£//*(* limn) sill k ipiMi sill t' t/> 0 , 

hi, ii 

I 11,0 


since : 


sin k ipim, sin f, HII | l*olynoin in oost/i H(H |, 


i.e. after summation, the terms with + n and - n cancel out. 
In addition: 

Ih il- Until) COSA 1 I '/I,,,,, COS A - m - o 

ih,ii 
i 0,0 

for k — 1,3,5,..., 

since cos (21 — 1) ^mn = cos '^mn /Polynom in sin^ '^mh/ 
i.e. terms with + m and -m cancel. 

We have : 

cos i (2 ,r) - - 2 l 1 cos* 2 * — 

i i li — :i\ 

- ( (i-:{)2' »cos' a2,«-(- 2 I ( lx 

X 2 (- 5 cos 4 '* 2x 1- i 


where : 


cos 2 x - cos 2 * a: - sin 2 a: (;« 2 - » 2 )/(»» 2 +» 2 ). 


If 1 is odd, i.e. 2i is not divisible by four, the terms cos L2x. 
reverse sign upon exchange of m and n, then the sum is equal to zero. 
If i is even, the signs are retained. So from equation 57: 
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and we immediately have, by coefficient comparison: 
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The coefficients from equation 6l and 62 are substituted into equation 


55 and we obtain: 
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Thus, the two scalar potentials are determined, except for one amplitude 
and propagation constant. 

7.3 Determination of the Vector Potential 


We set the emitting fraction of the vector potential equal to that 
of the individual fiber: 

A - ■(». <. 0 ) 

■ ax r: . (iso 


with 


This equation is axial symmetric. 

For the constant field the wave equations apply: 
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Since according to our assumption, the scattering field is 
Illustrated as a super-position of the emitting field, it follows 
that A'' 0 , since it cannot be excited by A ^ 

(By vector addition of all emitting terms, only components in the r, ip 
plane can exist and no z- component). 

As a general solution we select. 
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For the generated speeds of sound (only the fraction of the constant 
field) we have: 
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The field period is periodic in y- by tt/2: 

A"(<p -|- tc/2) - A" (<p). 

The field is symmetric at jh =0 (and because of the periodicity, also 
at <r • » n /2 ) . 

Therefore we must have : 


0 Mr 0) o 


and thus 


■I "tv 1 °) ■». 


and from this we obtain: 

l‘i 'h i i I A,- -h i 0 , 


from which follows: 

i>;; , /<;' • i>\’ ■- o f Qr i= 1 > 2 » 3 > • • • ( 74 ) 

2) Because of the periodicity, all terms where i is not divisible by 
four drop out . 


So as an equation meeting the conditions of symmetry and periodicity 
we have : 


A" ■■ 0 

ri ^sin 4 * 7' M u •/■!(* 1 (fr *') -t- 



*' ft 

+ 1 1 {<'>' '■) j t 

(75) 

yl” ~ c 

r * £ cos 4 i 7 ' [/I >h 1 - 1 (e V r) — 



— «4 , •/■! M I ('V>')l , 

(75) 

/i" - 0, 


(77) 


Equations 75 and 76 must still meet the divergence condition: 


div A" -- A r 

I / , 0 

A r + o r ' '' 


5 .H , * ® l"--0 

+ , V' 


or 


r(! or* ( t 'v , ‘) I* 

+ «4{'A(+1 (<■'>' r ) (’ 4~ ■**)] ' 

^ X' r . (< ( ’ 2 ^ sin 4 i 7 | /1 4 ; <•> '7 | i ■ i ( f i’ r ) 
Oc r f ,, , 

-I- i t • 


thus we have: 


( |ivyl"=o '' 2 ^ {sin 4 / 7 ‘| If (<A|,- -1 («»-»') 

I 

— ./.it 1 ( /■' r I')) I - I’ii <i.t f(’ 7 if 1 1 [t'r * ) i 

•I- * I •/■III I ('■>’■)) 1} • 
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because 


V .l,,(z) pi \( z )< 

1 z 

i- [ -i J p i(>) 


from this follows: 


div A" 


» n r v ^ Kin ‘l itpJn(fvf) X 
X ( — A,\j I- «',',) "0 


and from this we obtain: 


So as a solution we obtain the following expression which meets the 
condition of symmetry, periodicity and divergence: 


8 i A'^ sin -I i 7'. /.u ('■>'') , 

( 71 >) 

2J|,CO«-l 1*7 -1 ) - 

( 80 ) 

« 

( 81 ) 
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Now the constant field 79, 80 should originate as the superposition 
of the emitting field 64 of all individual fibers. These fibers are 
added vectorally. 

The emitting field of the fiber with the indices m, n generates 
the component A' at the plotted point P. This is split apart in the 
coordinate system (r, y )of the fiber with coordinates (0,0) into a 

A// A *' 

radial component and an azimuth component A v 

Prom Figure 2 we can see: 

« /? - -j2 , 

A" A' oosa A' cos {ft — it 12 ) A' sin/? , 

A” --- - A' sin a - - A' sin (/? • tc/2) A' con/? , 
sin/? ( /*’,„„/«) sin (»/»»«,* -7) (Kinossiit/.) , 
n- r- | li~ lw -- 2 r It m „ cos(iy>„ )n — 7O 
2((;osinnss;\ty,) . 

Key: 1-sine law 2-cosine law 


18 



Figure 2: The addition of vector potentials of the fibers (m, 

plotted point P. 

From this follows : 

cos// I' I sin'-// 

I 

w k />’ nni <•"*('/’«»« V )l • 

Thus we have : 


- /*’»«« «i" (ipiim v) 

' \- -|- l(~„„ •- 2 /' H m ,i cos (»/),„„ '/ ) 

(82) 

- /' — l{„i„vux(ii’i,ii, - <(‘) 

‘ 1 |/ -I- -• 2 >• cos (»/>„„, — 7 *) 

(82) 


thus : 


Now we should have: 


A 1 O ri lilllf x ( SJ ) 

X (/■> I V" + - 2r A“m»c"s(i/)„„, - 7 )) • 


,1 


/ / 
I* 



. 0.0 




IH . H 

1 0,0 



n ) at 


g] 



or for the radial component: 


I 


<■> r 


4 ,- si n (4 i f ) ./, K (f ,, /•) 

f 


(85) 


/4 '’ j >•! • 
} 9 

- 0.0 


/f i,i „ sin - <j ) 


^ mil 


-> H mn COS(y>,„„ — 7 ) 

■' "i ( Cv K "t I'mii 2 r //|h W i'os (»/>„,„ — i/ )) 


\ 


and for the azimuth component: 


i 

iii is 4 ifJ\i(f r r) ~ 

m 

v 

r limn <,(,s (V ’wh — '/ ) 

\ ^ 
~- f r) 

lj V /_, | 

in.)) 

| )- H~ m 2 r limn <ii)s(y wn 

I 0,0 

v//p (fl ,| 

, " 1 Ai",„ — 2 r ll t ini , ' ns (i/ l nui 

- r» • 


According to the Gegenbauer addition theorem, from the right side of 
equation 85 we have: 

9 // 00 

f|. /' * 0 

I 11,11 

x sin|(l I* k) iyu mi y )l ~ 

• 2/i ’" Sc I- *)•/..*(»>*•) 2* 

A- 0 

1 0.0 

X /^»«h) |sin(l +A-)v»»uiCOH(l -|-A:)y- 

~eos(l I A-) yi mN sin(l -|-A-)yj. 


Like the scalar potentials (see equation 59ff) the terms with 
sin(l + k) ^mn and terms with c-os(l + k) tf'mn, cancel where (1 + k) 
is not divisible by four. 

With (1 + k) = 4i from equation 87 we have: 

2 H'y 

„ „ 1 ^4 ((■,'*•) sill 4/y ^ II X 

-• (' fa, 

1 0,0 

( f 'r It III h) I'OS'I ? y,,,,, (KS) 
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And the coefficient comparison between equation 88 and left side 
of equation 85 Immediately yields: 

'hi ~ A'r //^ (<•> limn) ww(4 i y>inn) (NO) 

111,11 
t 0,0 

i— 1.2. :» 

Is not yet determined since the terms for 1*0 disappear on both sides. 


Prom the right side cf equation 86 according to the addition 
theorem we have: 


2A r V(l I 

r r* k o t n a 


i 0,0 


I'mii <'" s (y'Hi« <1 ) . ... . 

„ . , . «i»(* I £)('/'<«« '/)• 

OK)) 


/62 


A rearrangement which permits a direct coefficient comparison 
of r and 7 is not possible. Since A" is a term independent of angle, 
for A^ a coefficient comparison identical in r is sufficient for the 
angle independent fraction of equation 90 and the left side of equa- 
tion 86. For this, only the term for k = 0 or i = 0 is taken from 
equation 90 and the left side of equation 86. Then we obtain: 


2 1 C, 


2 'lo JiM “ 7 Me*r)'£ll'*\ ef ll mn )A 

* V f , 1 , „ 


III. II 

1 0.1) 


, . >' I'm 11 <'os(y>, H „ - 7) 

/*’ Mil 

2 K\ 

.. 'h (>!•>') V /q„„) 

wTh 


1 0.0 


k/tfww - cos i/», s 71 I sin i/i WH siiu; | . 


The terms with cos \b disappear in the sum due to alternation of n 

Mnn 

with -n. We then have: 


-■I., 


It Hill) 

' inTn ^ I'll III II 

1 11,11 


(!H) 
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From equations 64, 
potential A =. A* 4* 


79, 80, 81, 89 and 91 we obtain for the vector 
A": 


i/.S) 


- 2 '• 


/> <■ 


{•• /sin 4 / 7 ./,u(» •' 


1 ( 


(««) 


l‘nni) <'<>s •( / , 


W,M 
J 0.0 


<U(r t < h z) K[.v {//<-» (/•>/■) | 2 . /,(,>,') < 

2 2‘’<>*(4 m/)./, ,(/>/■) 

1 1 


^ #« o 

in.it * i' inn 

1 0.0 


^ H m I* mu) I * y'lmil 


Hl.lt 
t 0,0 


(»») 


Thus the vector potential is determined, except for one amplitude 
quantity . 

7.4 Comments 


For the constant field we have permitted an angle dependence, the 

emitting field of a fiber has been assumed to be axial symmetric. This 

is not entirely correct because the fields affect each other and an 
angle-dependent scattering field will always result in an angle-de- 
pendent emitting field. The solution presented can be considered to 
be an approximate solution. As a next step, we would have to proceed 
from an emitting field described by equation 63 , 92, 93 and to perform 
perturbation calculations with equations 79 and 80. This is not pos- 
sible either analytically or numerically. Within the framework of 

the achieved accuracy, we must assume that at least in the vicinity of 

the fibers an angle independent field exists. This is achieved by 
neglecting the last (angle dependent) sum in equations 63 and 93 and set- 
ting equation 92 equal to zero. This procedure is justified by con- 
sideration of the order of magnitude of the individual summands. 

O 

In equation 92 A^ disappears for r — » 0 on the order of 0(r ), 
likewise the last sum in equation 93, whereas the first sum in equa- 
tion 93 increases with 1/r. 


To formulate the boundary conditions at the fiber surface we 
therefore use the following descriptions of the potential: 




W,ll 
< 0,0 


v . //',■> (Cr «*») 
*'l' /'’mm 


Iftjf 
> 0.0 

A r A,. 0 , 


(! 0 ) 




The boundary conditions to be met are the same as for the Individual 
fiber. 

8 . Integration Method 

In the last section we Illustrated the scattering field as the 
discreet sum of the emitting fields. We expressly used the regular ar- 
rangement of fibers. But In a real absorber this regularity does not 
exist. We could account for this by Introducing an average fiber 
density and Illustrating the scattering field as an Integral over the 
emitting terms of the fibers distributed uniformly about an average. 

Let N be the number of fibers per surface unit. Then g is the 
amount for one fiber and NgdF Is the contribution of the surface element 
dF. 


In Figure 3 we determine the contribution of the source region /63 
about Q at the plotted point P. We again use equations 42 to 44 
for the individual fibers as the emitting field of the fibers: 

0 e ,oc : (^p, a 

A,, - /£,, //'f* (/•,,«)(> rc , 


m 
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Figure 3 : Integration relationships 

The surface element becomes 

(IF- FrtFdyj. 

From this results the contribution of the source region about Q at 
the plotted point P for the scalar potential: 

<10", - c r: N Fd F <ly> ■ 

X I I **" 2 V')) 

and from this : 

\ln oo 

K* ■ I' f« s ri NinCiWs 

0 <f 

X (r 0 , « |'V- h F» - 2 r A* « w (y« ~ f)) (IF. (»«) 

According to the addition theorem: 

H\\ H ,! o, a | | - F- - 2 /' F cos (yi ■— yi)) , 

•Ai ((V), a ( ) / /!r ) ((7., a F) ( 

00 

2 ^ >/« (i-’q, a f) F) cos|« (y> - yi) | , 


11 I 


Upon integration over \p the terms of the sum over n cancel because of 
the cosine. Therefore from equation 96 we have: 


n , ~TT 00 

0 t) ,,=(> rt J„(,- 0ta r) | dyi f F //},-> (,y. « F) (I F 

6 
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and from this : 


‘/V 




r t ,y v ,,, 

''o,oi ”« „ ,i -Al (/•(!, a >') . 




( 07 ) 


For the vector potential we must perform vector addition. Prom 
the discussion on radial or azimuth components in equations 82, 83 and 
84, we can immediately write the constant fraction: 


A' r ' (r,f,z) - 


7 I n 


// 2 siii(yi ■■-</■) 


.. t ' I* .» #1 ('III V I 

o (<I V > , „ , .. . . x 

„ ,/ | *•- I K* - 2 >' /> «>s (f - y) 

y //',-*(/> I 7 *'* I /i a — 2r /fcos(i/> - </ )) dll , (08) 

'I T< s) * 

, i i n f >• — // cos (v» — <r) 

o A,, iV /l> | t r2 l{ > _ 2 rlt cos ( y* - V-) * 

/ I'/ " f />’“ ~ 2 r // cos ( y> - y)) <1 /( . (00) 


The integration limits with regard to if/ are variable as long as 
integration proceeds over a full circle. The selected shape permits a 
simple transformation: a y> - y . Thus from equations 98 and 99 

we have : 


„ „ /i 2 sin a 

X //'-> (o |/r 2 + « 2 - 2 r Jl cos a) (1 It , ( ( 00) 

r — Jl cosa 

X // ( f ) (s ,, I'V-’ -I- 11* - 2 »• R cos a) (1/i . (101) 


In equation 100 the integrand with respect to ct is an odd function, i.e., 

A" - 0. 
r 
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In equation 101 the integrand with respect to a is an even function. 


According to the addition theorem we have: 

| //•“’ -|- It- — 2r It i!os «) - 

" .. /. ^ I- Hi\k(tfll) (V 2 I- K 2 — 2 r It cosy. (102) 

t i > "i o sina 

And thus from equation 101: 


J. f ■ I’l It]' K oo n 

A", " f ((>• — It cos a) X (100) 
ii 

». sin A*a , , 

X yk'J k ^,i)ll[r ) (F f lt) . <latl It. 

sina 

In equation 103, If r| I, , then j ./*(/■>,•)) very quickly goes 

to zero with Increasing k. In addition, the term 'sin^a/Kina, 
oscillates with increasing k more and more so that upon integration, 
only a slight contribution is provided. Of the sum appearing in the 
integrands we therefore consider only the term with k= 1. We then 
have : 

„ 4 v n /<]',, N °?p 

A, t ■— | f {r — It cosa) X 

f V r ,) 0 

x ./l (/ ,./•) ('V lt)Aou\li . 


Upon integration over a the first term of r - R cos a yields 
the factor 7 r, the second term disappears. It becomes: 

.l” A-*., I'zh/ A r * ( ^ »’ ^ ) 


A, 4*, 

if ,v 


, //!r*(ivrf) 

4 7r l: li,.Ntl- lv , -I i (r,r) 
V'r"r 


Thus we obtain a description of the potentials 


<1*0, a ( r > z ) l' j o,y. v 1 ^1 (t'e, a 0 — 

// ( , U, (/'V> a rf) 

- 2 nN (I- J 0 {l> 0 ,*r)\, 

K(l, * « 

/l^(r.s) ' > |// ( ,- ) (K ) ,r) f 

, , , v //.fffrd) , , .. 

■I- - 4 nN(f- J\ (/••,,/•)! , 


(<><l) 2 

A r = A z := 0 . 


With these descriptions the boundary conditions at the fiber surface 




1 


I 

are again fulfilled. They correspond to equations 94 and 95 for the 
model with quadratic, regular fiber arrangement and are comparable 
with these equations. 


We have now described the fields of the compressional , temperature 
and viscosity waves in the environment of a fiber. Both in the complete 
field equations 6 3 , 92, 93 as well as in the simplified equations 94, 

95 for a regular quadratic fiber arrangement and in equations 105 , 106 
for a homogeneous fiber distribution, the amplitudes Ep <eitl; • and the 
axial propagation constant P are unknown. These are, as already shown 
in Section 6 for the individual fibers, determined through the boundary 
conditions of equations 33» 34, 35 at the surface of the reference 
fiber (m, n)= (0, 0), i.e., at r~ a. 


Prom a comparison of the equations 105 and 106 by the integration 
method with the results of the summation for regular fiber arrangement 
we now know that the component A” (r , y. , z) of the scattering field 
quite disappears. This provides an additional argument for the simplifi- 
cation performed in Section 7.4, mainly, instead of using equation 92, 
to simply set A^Cr,'/ 1 , z) = 0. We also find a relation between equa- 
tions 94, 95 and equations 105, 106: 

t 0,0 


r(ty.arf) 

*'t>, a <1 


“ If It lilt) 


2 7t JV <l~ 


tv," 
i 0,0 


(< r '/)" 


In the complete field equations 63 , 92, 93 we see in the last terms 
the influence of the spacial arrangement of neighboring fibers on the 
total field in the vicinity of a reference fiber. Prom the derivation 
of these terms as well as by its form, the influence of the order of 
symmetry is seen. If this increases, then in accordance with the factor 
in front of the summation index i, i.e., the order, of the first cylinder 
function increases in the terms and the interval between the orders 
becomes larger . The sums in the terms converge more quickly. 
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At this point we can account for the microstructure of the real 
absorber by introducing the statistically most frequent order of 
symmetry and by performing the summation of scattering fields of 
neighboring fibers only up to a certain interval and by applying the 
integration method to more distant fibers. 

In part 2 of this work we will utilize the randomness of symmetry 
selection to derive a more simplified model. 

10. Conclusions 


One objective of absorber theory is the obtainment of a calcu- 
lation method for the propagation constants and for the wave resistance 
of the absorber. Another objective is to gain knowledge about the 
influences of the individual structural parameters and material data 
on these values. Both objectives will be pursued in part 2. 

Thanks are due the German Research Society at this point for 
financial support for the work. 
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